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1. INTRODUCTION 


a. Motivation for .Iny.es.tlaatiog 

Atmospheric motion is characterized by sizes ranging from 
planetary waves down to small scale turbulent eddies. Scale 
interactions are observed to transfer energy both up and down this 
meteorological spectrum. Understanding the complex interactions 
between scales remains one of the most challenging topics in 
meteorology. Mesoscale phenomena such as severe thunderstorms are 
influenced by both larger and smaller scales of motion. In their 
developing and advanced stages, thunderstorms apparently modify 
the wind and thermal structure of their surrounding environments 
through scale interaction or "feedback" mechanisms which presently 
are not well understood. 

Inadequate data for describing mesoscale weather phenomena 
have hindered research into scale interactions. The resolution of 
the routine upper air network is too coarse in both space and time 
to properly capture mesoscale weather events since observations 
are taken every 12 h at an average station spacing of 400 km. 
However, special mesonetworks developed for the Atmospheric Varia- 
bility Experiment - Severe Environmental Storms and Mesoscale 
Experiment (AVE-SESAME) (Alberty sL. , 7979 ; H i 1 1 3L 2k. , 1979 ; 
Barnes, 1981 ) are allowing researchers to probe, measure, and 
observe the mesoscale storm environment at a much finer resolution 
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than is routinely possible. Many investigations utilizing the 
special data are now underway to learn more about the complex 
structure of mesoscale weather phenomena and their interactions 
with the i«*»*»er scale environment. Our ability to forecast weather 
events will improve, for example, as our knowledge of scale 
interactions increases and the findings are incorporated into the 
numerical prediction models. 

One way to increase our understanding of thunderstorm- 
environment interactions is through an energy analysis of the 
storm's surroundings. The research described here concerns the 
environmental kinetic energy budget of an area of intense thunder- 
storms that occurred during the fifth AVE-SESAME period (20-21 
1979 ). By using meso g-scale data with an average station spacing 
of 75 km at either 1.5 or 3 h intervals, a detailed description of 
variations near the storms over central Oklahoma can be obtained. 
The energy processes are examined prior to, during, and after peak 
storm activity in order to better understand energy fluctuations 
related to the convection. 

b. Past studies 

Diagnostic studies using conventional rawinsonde data suggest 
that organized convective activity alters its larger scale 
environment. Numerical simulations by Aubert ( 1957 ) showed that 
latent heat release associated with convective activity produced 
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increased low level convergence and enhanced upper level diver- 
gence. In addition, height rises were found to occur in the upper 
troposphere, whereas height falls occurred near the surface. 
Using Conventional rawinsonde data and satellite imagery, Ninomiya 
(1971) detected a mid- tropospheric warm core, strong low level 
convergence, and increased upper level divergence in the storm 
surroundings, A mid- tropospheric jet maximum formed northwest of 
the convective area in response to a strong thermal gradient. As 
the storms developed, a low level jet intensified nearby. 

More recently, Maddox ( 1980), Maddox e.t al. ( 1980) , and 
Fritsch and Maddox (1980, 1981a, b) used enhanced, infrared satel- 
lite imagery and conventional rawinsonde data to define and 
describe the Mesoscale Convective Complex (MCC). MCC's are large, 
continuous, nearly circular areas of intense storms that last for 
at least 6 h. The size of the MCC cloud shield may exceed that of 
an individual thunderstorm by more than two orders of magnitude, 
They are thought to modify their larger scale environment to the 
extent that the evolution of future synoptic- scale systems is 
affected (Maddox, 1980). Environmental conditions near the com- 
plexes include areas of mid- tropospheric warming and ascent, 
enhanced low level convergence and upper level divergence, regions 
of upper level ridging and anticyclonic flow, and the formation of 
jet maxima on their poleward sides. Numerical simulations by 
Fritsch and Maddox (1981b) have duplicated some of the larger 
scale environmental modifications associated with areas of convec- 
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tive activity. 


In addition to the synoptic- scale studies, mesoscale investi- 
gations of storm-environment interactions include those by Fan- 
khauser (1969, 1971, 1974), Sanders and Paine (1975), Ogura and 
Chen ( 1977 ) » Sanders and Emanuel ( 1977 ), Ogura and Liou (1980), 
Moore and Fuelberg (1981), Vincent ji. (1981), Anthes jjl,. 
(1982), Park and Sikdar (1982), and Vincent and Carney ( 1982). 
These diagnostic studies have confirmed the findings that were 
observed at the synoptic scale. Although, it is believed that the 
thunderstorms produced many of the observed environmental varia- 
tions, mote detailed investigations are needed to understand the 
mechanisms involved. 

The literature contains relatively few studies describing the 
energetics of the convective storm environment. Most previous 
diagnostic kinetic energy investigations have focused on 
synoptic- scale systems such as extratropical cyclones. Those stu- 
dies that have considered the energetics of the convective 
environment suggest that large areas of storms modify the sur- 
rounding synoptic- scale energy balance. Vincent and Schlatter 
(1979) proposed that deep cumulus convection could supply kinetic 
energy to the synoptic- scale flow. In their study involving the 
remains of tropical storm Candy, areas of kinetic energy increase 
and positive dissipation agreed closely with locations of the con- 
vection. Robertson and Smith (1980) investigated the energetics of 
two extratropical cyclones that produced the Jumbo (3-4 April 
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1 97-4 ) and Palm Sunday (10-11 April 1965) tornado outbreaks. Using 
synoptic- scale rawinsonde observations every 12 h, they found 
major energy activity associated with the cyclones that produced 
the outbreaks. 

With the completion of the Atmospheric Variability Experi- 
ments (AVE) and the AVE-SESAME projects, studies incorporating 
higher resolution data now are being reported. Fuelberg and Scog- 
gins (1978) used AVE IV rawinsonde data having a synoptic- scale 
spacing but a subsynoptic- scale time interval to compute the 
kinetic energy balance of the Mesoscale Convective Complex 
described by Maddox et al. O98O). Rather tranquil conditions 
existed for the analysis region as a whole, whereas very active 
energy processes occurred within the areas of convection. The 
storm environments were characterized by generation of kinetic 
energy by cross-contour flow, horizontal flux divergence near jet 
stream level, upward transport by large scale ascending motion, 
and transfer of energy from resolvable to non-resolvable scales of 
motion (negative dissipation). They also noted that maximum 

energy conversion and transport occurred near the times and loca- 
tions of greatest storm activity, thereby suggesting that energy 
variations were due, in part at least, to the convection. 

Using National Weather Service (NWS) rawinsonde observations, 
Fuelberg ct al . (1980) examined the synoptic- scale kinetic energy 
balance during the deadly Red River Valley tornado outbreak which 
coincided with the first AVEASESAME period (10-11 April 1979). 
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The effects of upper level jet intrusion dominated the energy 
budget for the entire analysis region. Horizontal flux conver- 
gence and cross-contour destruction of kinetic energy associated 
with supergradient flow were the major processes. In contrast, 
the energetics of the storm environment included generation of 
kinetic energy and positive dissipation. Fuelberg ct al . (1980) 
could not isolate possible feedback mechanisms from the storms due 
to strong, pre-existing dynamic forcing. They hypothesized that 
subsynoptic- scale data would be useful for resolving storm- 
environment interactions in these situations. 

A subsynoptic- scale kinetic energy study of the Red River 
Valley tornado outbreak was recently conducted by Fuelberg and 
Jedlovec (1982). The combination of NWS and special rawinsonde 
observations provided a superior resolution of smaller scale 
features than was obtained by Fuelberg et al. (1980). Once 
again, it was found that the energetics of the storm environment 
differed greatly from that of the overall analysis region. Magni- 
tudes of energy processes were larger at the finer scale of reso- 
lution, and effects of the convection were more apparent. 

The only kinetic energy budget studies of the meso g-scale 
storm environment appear to be those by Mclnnis and Kung ( 1972), 
Kung and Tsui (1975), and Tsui and Kung (1977) who utilized rawin- 
sonde data from the National Severe Storms Laboratory (NSSL). 
With an average station spacing of 85 km and observations at 1.5 - 
3 h intervals, the data provided a detailed description of energy 
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variations within the mesoscale storm environment, The convective 
environment was energetically active and characterized by intense 
generation and dissipation of kinetic energy in the upper layers 
of the atmosphere, Ch the other hand, periods without storms con- 
tained strong destruction of kinetic energy by cross-contour flow 
and transfer of energy from subgrid to grid scales of motion. 
Magnitudes of meso 8- scale energy processes varied widely depend- 
ing on the presence and strength of the convective storms, but in 
some cases, values were an order of magnitude larger than 
corresponding processes on the synoptic scale. 

c. Objectives 

The objective of this research is to investigate the kinetic 
energy balance of a mesoscale storm environment as a way to 
increase our understanding of how areas of convection modify their 
larger scale surroundings. Special emphasis is placed on the 
energy balance in order to explain observed wind variations. Meso 
3-soale rawinsonde data from the fifth AVE-SESAME period are 
used. In contrast to mesoscale studies by Kung which emphasized a 
collection of days, the current investigation is a detailed 
analysis of a single 24 h period. Horizontal maps of energy 
parameters and other diagnostic results that have not been 
presented previously at this scale will be described. 

The specific topics that are investigated include 
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1 ) descriptions of wind and kinetic energy variations within the 
analysis region, 

2 ) relationships between temporal variations in the kinetic 
energy budget, wind variations, and the convective activity, 

3) agreements between horizontal patterns of energy parameters 
and locations of the convection, 

4) a comparison of the energetics of the same analysis region 
derived from different scales of resolution, and finally, 

5) an error analysis on the derived energy parameters. 
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2. METHODOLOGY 


a. Theoretical develop gent 


Both synoptic- scale motions and subsynoptic- scale motions 
outside of convective cells are quasi-horizontal and virtually 
hydrostatic. Therefore, the theory used in large scale energy 
analyses can be extended to meso 8-scale studies by utilizing 
higher resolution data (Kung and Tsui, 1975 ). This study uses the 
kinetic energy budget equation for a fixed, limited volume in iso- 
baric coordinates which is given by Ward and Smith ( 1 976 ) as: 

|| = //-V.V.J. - //V-kV - + //V-F + fk (1) 

ct dp 03 t , 

(A) (B) (C) (D) (E) (F) 

where 

If - fff dx dy dp. 

If is the horizontal wind velocity, 

w is the vertical motion in isobaric coordinates, 

2 2 

k = (u +v )/2 is the horizontal kinetic energy per unit mass, 

K = Ilk, 

0 s gz is geopotential height, 

f is the frictional force, 

A is the computational area, and 

is a subscript denoting surface values. 
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Local changes of kinetic energy , term [A) above, are due to 
five processes. Term .(B) represents generation of kinetic energy 
due to cross-contour flow over geopotential height gradients. 
Terms (C) and (D) denote horizontal and vertical components of 
flux divergence, respectively. Term (F) represents variations in 
energy due to changes in mass of the volume being studied. Its 
value was found to be several orders of magnitude smaller than 
those of the other terms; and, therefore, it will not be con- 
sidered further. Finally, term (E) conceptually represents ther- 
modynamical and mechanical frictional processes. This term cannot 
be treated explicitly and is therefore computed as a residual to 
balance the other terms of the equation. In this way, the resi- 
dual term not only represents frictional processes, but also 
transfers of energy between resolvable and unresolvable scales of 
motion, together with possible error from the calculation of the 
other terms. It is commonly called the dissipation term. Since 
the resolvable scales of motion are known, this term allows one to 
infer the importance of motions which are not readily detected 
with the input data and computational procedures being used. 

h Data 

The case being studied is the fifth day (20-21 May) of AVE- 
SESAME 1979. Rawinsonde data were gathered from twenty sites in 
Oklahoma and Texas at either 1.5 or 3 h intervals, i.e., on 20 May 
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at 1100, 1400, 1700, 2000, 2130, and 2300 GMT and on 21 My at 
0200, 0500, 0800, and 1100 GMT. The average station spacing was 
approximately 75 km (meso 3-scale.) which is similar to that of 
data sets collected by the National Severe Storms Laboratory 
(NSSL) ■ In addition twenty-three National Weather Service (NWS) 
sites surrounding the mesonetwork also collected data every 3 h. 
These stations had the conventional synoptic- scale spacing of 400 
km. Figure 1 shows the location of the meso 3 (storm scale) net- 
work in relation to the routine rawinsonde stations together with 
a detailed view of the special sites. Soundings were terminated 
at 100 mb, and no data were gathered for Norman, OK (0UN). Mimes 
of the special sites are given in Table 1 . Reduction procedures 
used to process the rawinsonde data are described by Fuelberg 
(1974), while tabulated results at 25 irb intervals are given by 
Sienkiewicz et al . (1981). 

It must be emphasized that we are investigating energy fluc- 
tuations of the gtorm environment rather than individual convec- 
tive cells. Although the special mesonetworks of AVE-SESAME 1979 
have a much higher density than ordinary upper air sites, even the 
storm scale network still is too coarse to resolve details of par- 
ticular cells. According to Sienkiewicz (198I), only four 

sondes from the special network exhibited large variations in 
ascent rate which might indicate penetration into active convec- 
tive cells. A careful time and space continuity check of the four 
soundings, however, suggested that they also were representative 
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Fig. 1. Rawinsonde networks participating in AVE-SESAME V. The 
top diagram shows the synoptic— scale network while the 
bottom gives the meso g-scale network. The storm scale 
domain is superimposed on the large scale domain in the 
top figure. Outlines (bottom) indicate the 15 X 13 
analysis region (outer), the overall energy averaging 
region (middle), and the subvolume encompassing the 
upper level wind couplet (inner) . 
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Table 1. Names of the special storm scale radiosonde 
stations participating in AVE-SESAME 1979. 


ADA 

Ada , OK 

LTS 

Altus, OK 

CAN 

Canadian, TX 

CHE 

Cheyenne, OK 

CHK 

Chickasha, OK 

CDS 

Childress, TX 

CSM 

Clinton Sherman, OK 

EMC 

Elmore City, OK 

FSI 

Ft. Sill, OK 

GAG 

Gage, OK 

HEA 

Healdton, CK 

HEN 

Hennessey, OK 

HNT 

Hinton, OK 

TVY 

KTW, OK 

MTV 

Mountain View, OK 

OUN 

Norman, OK 

SEL 

Seiling, OK 

SHM 

Shamrock, TX 

SUD 

Stroud, OK 

SPS 

Wichita Falls, TX 
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of the storm environment, i.e. f no unusual characteristics were 
noted. Moreover, the objective analysis procedures described next 
should greatly filter any non- environmental effects that escaped 
detection. Thu°. data from the four soundings could be used in 
the later computations. 

As an aid in the diagnostic study, radar summary data from 
the National Meteorological Center (NMC) and GOES satellite 
imagery obtained from NASA's Marshall Space Flight Center (MSFC) 
also were utilized. 

c. Data handling procedures 

The use of subsynoptic- scale data, may violate several 
assumptions that normally are made as part of large scale data 
handling procedures. Furthermore, since subsynoptic- scale resolu- 
tion depends on having as much data as possible, it is desirable 
to generate values for missing soundings when feasible. Conse- 
quently, several modifications were made to the usual synoptic- 
scale techniques. 

Rawinsonde data collected during AVE-SESAME V were generally 
found to be of good quality. Nevertheless, several problems 
remained in the final tabulated product. First, data from ten 
storm scale soundings appeared questionable because their values 
did not show reasonable time and space continuity (Table 2). 
Therefore, these soundings were removed from the data set. 
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Table 2. Mesonetwork soundings removed from the 

data base because of abnormal characteristics . 


Station 

Date/Time 

Cheyenne, OK 

20 May/ 1100 GMT 
20 May/ 1400 GMT 
20 May/ 1700 GMT 
20 May/ 2000 GMT 
20 May /2 130 GMT 

20 May/2300 GMT 

21 May/0800 GMT 

Seiling, OK 

20 May/1700 GMT 

20 May/2130 GMT 

21 May/ 1100 GMT 
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Secondly, some releases were either missing altogether (11 
occurrences), terminated early, or contained data gaps at various 
levels because of sonde failure or other difficulties. In addi- 
tion, although most of the sondes were released within 15 minutes 
of each other, a few releases were delayed by as much as 30 to 90 
minutes. Thus, non- simultaneous releases combined with variations 
in ascent rates caused the sondes to reach given pressure levels 
at different times. Fankhauser (1969) has noted that for mesos- 
cale studies, the data must be adjusted so they describe a partic- 
ular level at the same time. 

To deal with these problems, analytical procedures similar to 
those of Fuelberg and Jedlovec (1982) were utilized. To account 
for non-simultaneous sonde releases and variable ascent rates, 
data at each level were adjusted to a common time using a linear 
interpolation scheme similar to that of Fankhauser ( 1969). Since 
the scheduled release times were 55 minutes prior to the standard 
hour (e.g. f 1105 for 1200 GMT), all data were adjusted to the 
prior hour (e.g., 1100 GMT) and will be denoted as 1100 GMT, 1400 
GMT, etc. soundings. To generate values for missing data, linear 
time interpolation at individual stations was used if two condi- 
tions were satisfied. First, interpolation was not employed over 
a time interval greater than 6 h, and second, extrapolation was 
avoided by not using the procedure at the first or last times. 
Approximately 15% of the total soundings used during the 24 h 
period were created in this way. Table 3 gives a list of avail- 
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able data before and after the interpolation technique. As a check 
on the final data set, time series and constant pressure plots 
were computer generated and hand analyzed. Values that did not 
appear valid were removed and not used in any further calcula- 
tions. 


After obtaining the time consistent data set, the storm scale 
values were objectively analyzed onto a 15 x 13 grid using the 
Barnes (1964) procedure. A grid length of 25 km and a scan radius 
of 150 km (2 times the average station spacing or 6 grid lengths) 
were used. Four iterations of the scheme produced the desired 
resolution. Gridded fields of height, temperature, winds, and 
dewpoint were obtained at the surface and at 50 rrb intervals from 
900 to 150 mb. Too much data were missing at 100 mb to permit 
analysis. Since average sonde drift at 150 irb was approximately 45 
km, nearly two grid units, sonde locations at individual levels 
were determined as part of the analysis procedure and used instead 
of surface locations, Winds at individual 50 rrb levels were aver- 
ages at that level and 25 irb above and below in order to suppress 
inherent data errors. A mild filter was applied to all gridded 
fields (Shuman, 1957) to diminish small waves that could not be 
adequately resolved. 

With an average station spacing of 75 km, unambiguous resolu- 
tion is restricted to features with wavelengths greater than 150 
km (two times the station spacing). Therefore, phenomena with 
smaller wavelengths should be greatly suppressed. In addition, 
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since the analysis region is only about 400 km on a side, larger 
features are only partially resolved by the special network. Fig- 
ure 2 shows the response curve for the storm scale study. The 
final grids retain approximately 32$ of originally resolved ampli- 
tudes at 150 km wavelengths, and 82% at 225 km wavelengths (two 
and three times the average station spacing, respectively). 

The Barnes ( 1964) analysis technique also was used to prepare 
gridded fields of the synoptic- scale data. Using a mesh length of 
127 km, a scan radius of 371 km, and four iterations, the values 
were objectively analyzed onto a 15 x 13 grid. Procedures were 
similar to those described previously for the special network. 
The resolution profile for the synoptic- scale data also is shown 
i n Fig. 2. 

d. Commitational procedures 

Since details of the procedures used here are given by Fuel- 
berg and Jedlovec (1982), they will only be summarized in this 
report. The kinematic method was used to determine vertical 
motion (oj). The lower boundary condition, normally assumed zero 
on the synoptic scale, was computed as the sum of three processes: 
local change of surface pressure, cross- isobaric flow, and 
terrain- induced vertical motion. A limitation of the kinematic 
method is that inaccuracies in the winds accumulate during the 
integration of divergence so that unrealistic magnitudes of u 

19 


J 




Fig. 2. Response curves for the meso 3-scale (dashed) and 

synoptic-scale (solid) gridded fields. Wavelengths 
are in km. 
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occur near 150 mb. To circumvent this problem, values of both to 
and horizontal divergence were modified using O'Brien’s (1970) 
adjustment scheme. This procedure forces the vertical motion to a 
predetermined value at 150 mb. Then, values of divergence are 
altered by amounts that would produce the modified to profiles. 
Our original goal was to set kinematic motions at. 150 irb equal to 
values obtained from the adiabatic method. Adiabatic motions can 
be obtained independently of the kinematic results and were 
hypothesized to be more realistic than an assumption of zero which 
is normally used. Unfortunately however, trial adjustments to 
adiabatic values did not produce reasonable results, apparently 
because the intense convection produced significant latent heat 
release that greatly violated the adiabatic assumption. Therefore, 
vertical motions were adjusted to zero at 150 mb in the usual 
manner . 

Terms of the kinetic energy budget equation were computed 
using horizontal and vertical centered finite differences at the 
17 levels. These values then were integrated over 50 mb layers by 
means of the trapezoidal rule. Dissipation was calculated as a 
residpal from the other terms in (1). The use of non-centered time 
differences for calculating, 3K/3t was investigated. However, 
this method is somewhat unsatisfactory since a complete energy 
budget is not produced at every observation time. Results showed 
that values of both 3K/3t and dissipation generally were not 
affected greatly by the type of time differencing employed ■ 
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Therefore, centered finite differences were used. The only excep- 
tions were at the first and last times and when the 2130 GMT ob- 
servation required an uncentered procedure. 

e. Error analysis 

Uncertainties in rawinsonde data are mainly due to random 
error, systematic error, and aliasing. The data handling and com- 
putational schemes described earlier should greatly reduce the 
effects of these errors on the energy budget calculations; how- 
ever. they will not be completely suppressed. Therefore, it is 
imperative to determine quantitative confidence limits for the 
derived parameters. Unfortunately, there is no effective pro- 
cedure to test for systematic errors; therefore, this type of 
uncertainty was not considered. In addition, data limitations and 
computational inadequacies such as truncation error were not 
evaluated. Ch the other hand, the effects of random errors on the 
energy analyses were documented using randomly perturbed data. The 
procedures utilized in this analysis, as well as detailed results, 
are presented in the Appendix. Briefly stated, however, the find- 
ings indicate that errors in rawinsonde data should not signifi- 
cantly affect the interpretations of results that follow. 
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3. WEATHER CONDITIONS 


Although the AVE-SESAME I period has been the subject of many 
studies because it coincided with the Red River Valley tornado 
outbreak, the fifth SESAME case (20-21 May 1979) also was a period 
of significant convective activity. During the 24 h period begin- 
ning at 1200 GMT 20 May, a total of seven tornadoes, one funnel 
cloud, and several hook echoes were documented within the overall 
synoptic- scale network (July and Turner, 1981). In addition, 
there were numerous occurrences of hail. Although no tornadoes 
were reported within the meso B-seale storm network, thunderstorm 
activity was intense, and rainfall amounts exceeded 3 in. at 
several sites. This section highlights weather conditions of the 
storm scale network during AVE-SESAME V; however, a more detailed 
account is given by Fuelberg and Printy ( 1 983 ) . 

Synoptic conditions at the beginning of the experiment (1200 
GMT 20 May) are shown in Fig. 3. At the surface, a relatively 
strong cold front extended from the Great Lakes into the Ohio Val- 
ley. From there it became stationary, winding through central 
Oklahoma and into New Mexico. Conditions in the upper troposphere 
were dominated by a cut-off low centered just south of Arizona and 
its major trough extending northeastward. A short wave trough was 
situated in the upper Midwest. Jet maxima were located over the 
Iowa-Minnesota border and south of the cut-off low over northern 
Mexico. Winds were generally light over Oklahoma and northern 
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Texas. Precipitation was occurring north of the front, covering a 
broad area of Missouri, northern Arkansas, and eastern sections of 
Kansas and Oklahoma. Echo tops of thunderstorms reached 17.1 km 
(56,000 ft), and rainfall totals of 1-2 ■*h» occurred in these 
areas during the early morning hours. The storms had moved east- 
ward, out of Oklahoma by 1500 GMT when they were situated over 
Arkansas and southern Missouri. 

The front gradually moved southeastward during the day, and 
by 1700 GMT, it extended from near Stroud (SUD) in northeastern 
Oklahoma to near Shamrock (SHM) in the Texas panhandle. The atmo- 
sphere south Of the front was quite unstable by midday. At 1700 
GMT, values of the Total Totals and K indices were generally in 
the lower 60’s and upper 30's, respectively, while Lifted Indices 
were near -6. This instability coupled with afternoon heating and 
frontal lifting provided the necessary conditions for intense con- 
vective activity. 

Radar summaries for the peak storm period are given in Fig. 
4. In addition, constant pressure maps for the meso £ -scale 
region are presented in Figs. 5-6 for 850 and 200 mb, respec- 
tively. Just prior to convective development (1700 GMT), the 850 
mb chart depicts a low over the Texas panhandle and a slight ridge 
over central Oklahoma. Winds were weak but exhibited considerable 
cyclonic turning. The 200 mb height contours were generally 
smooth, showing no major perturbations. Southwesterly winds con- 
tained very little diffluence. 
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Fig. 5. 850 mb height analyses in 5 m intervals (solid ) ; 00 

means 1500 m. Winds are in m s -1 with short barbs 
denoting 5 m s~^ . Isotherms (dashed) are in °C. 
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Convection first formed along the front near 1700 GMT, and by 
1935 GMT, a well defined area of storms extended from near Chil- 
dress, TX (CDS) to Ada, OK (ADA) (Fig, 4), The strongest cells 
were located between Childress and Wichita Falls, TX (SPS). Max- 
imum echo tops in this region reached 17.7 km (58,000 ft), and 
hail was reported near Ada. Figure 7 shows the extent of cloud 
cover over the mesonetwork at 2026 GMT, The intense convection 
appeared to be most widespread near 2135 GMT (Fig. 4). Radar 
reports from Oklahoma City (OKC) at 2330 GMT indicated that level 
4 thunderstorms (TRW++) still covered most of Oklahoma and north- 
central Texas. By 0135 GMT, light rain persisted over Oklahoma as 
the heavy storm cells moved eastward into Arkansas. However, 
other strong storms were nearby, situated just south of the spe- 
cial network along the Red River Valley. 

Significant changes in the atmosphere occurred during peak 
storm activity, and these are clearly seen on the constant pres- 
sure maps (Figs. 5-6). Surface pressures (not shown) fell over 
most of the network between 1700-2130 GMT. Greatest decreases of 
3-5 mb were located in the southeastern half of the data region. 
Consequently, height falls of 20-40 m at 850 nb (Fig. 5) are com- 
mon over southeastern Oklahoma during the 4.5 h period. This 
results in a pronounced trough from near KTVY (TVY) to near Chil- 
dress (CDS). Winds at several stations in the southeastern part 
of the network have shifted to a more southerly direction, thereby 
supplying warm, moist air into the convective area. At 200 mb 
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(Fig, 6), significant warming and ridging has occurred in close 
proximity to the storm area. Winds now show pronounced diffluence 
and strongly ageostrophic flow. Maps at 0200 GMT suggest that the 
region was returning to pre-storm conditions. Temperature and 
height patterns at 850 rrb still are rather disturbed; however, at 
200 mb, the ridging has greatly diminished, and diffluence is much 
less pronounced. 

Strong variability within the meso B-scale storm environment 
can be seen in pressure-time cross sections of horizontal diver- 
gence and vertical motion (Figs. 8 and 9, respectively). These 
diagrams represent area-averaged values for the innermost 13 x 11 
grid points of the total analysis region. The time series for 
horizontal divergence depicts the rapid development of pronounced 
low level convergence and upper level divergence after the convec- 
tion began shortly after 1700 GMT. Greatest values coincide with 
the times of peak storm activity. These findings are consistent 
with those observed on the constant pressure maps, i.e., strong 
inflow at 850 rrb and diffluence at 200 mb. Values decreased after 
0200 GMT as the intense storms were replaced by rainshowers. 
Strong mesoscale ascent developed in response to the wind flow 

(Fig. 9). Although maximum area-averaged vertical motion was -22 

-1 

pb s near 400 nb at 2130 GMT, horizontal cross sections for that 
time (not shown) indicate that values at individual grid points 
exceeded -60 yb s”"' ■ 

During the last half of the experiment period, the cut-off 
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Fig. 8. Pressure— time cross section of horizontal 

divergence (10 - ^ s“*) for the inner 13 x 11 
grid area of the total analysis region. 
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li 14 17 20 23 02 05 08 11 

Time (GMT,) 

Fig. 9. Pressure-time cross section of vertical motion 
(yb s -1 ) for the inner 13 x 11 grid area of the 
total analysis region. 




low in the upper levels moved slightly northeastward toward the 
Arizona-New Mexico border. The front advanced southeastward with 
continuing isolated severe convective activity nearby. By 1200 GMT 
21 May, the end of the period, the cold front had become station- 
ary over central Texas. Although storm intensities had decreased, 
convective activity continued over Texas, Oklahoma, Arkansas., and 
southern Kansas. 
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4. RESULTS 


This section will describe the kinetic energy balance of the 
meso 3-scale' network. Area-averaged energetics will be discussed 
first, followed by a description of its time variability. Next, 
horizontal maps of kinetic energy budget terms will be presented 
to explain flow patterns within the storm scale environment, 
Finally, the energetics of a subregion over the southeastern por- 
tion of the network will be isolated to diagnose possible causes 
of upper level wind variations that are observed during the times 
of intense storm activity. 

a. &ea-time averaged en emetics 

Terms of the kinetic energy budget equation were evaluated 
for each observation time at individual grid points. Integrated 
grid point values were averaged over the energy analysis region in 
Fig, 1. Area-averaged results for the ten times then were com- 
bined to produce the composite area-time averaged energy budget 
which is given in Table 4. 

Kinetic energy content for the surface to 150 irb column is 
11,99 x 10 J m t with the major contribution coming from near 
jet stream level. All levels experience a net decrease of kinetic 
energy during the 24 h period, with the vertical total being -5.41 

mmO 

W m « This local change is due to various sources, sinks, and 
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Vertical 

Total 11 99 -T 41 83 57 -G7 39 0 03 -41 6Z -1C 19 -33 20 



transports. Generation of kinetic energy by cross-contour flow, 

_ p 

which produces a vertical total of 83.57 W m , is the dominant 
source to the region. Energy generation occurs within most layers 
but is greatest above 400 mb. A secondary maximum exists I.. the 
planetary boundary layer. 

Horizontal flux divergence and dissipation are major sinks of 

-2 

kinetic energy. Vertical totals for these terms are -47.39 W m 
and -41.62 W m , respectively. Again, the processes are largest 
in the upper levels of the atmosphere. The vertical profile of 
dissipation is similar to that of generation, i.e., a maximum near 
jet stream level, weak dissipation in the mid-troposphere, and 
another maximum near the surface. Significant horizontal outflow 
exists in the upper troposphere (negative values), while weaker 
horizontal inflow occurs in the middle layers. 

To further investigate the physical mechanisms creating hor- 
izontal flux divergence, the term was separated into two com- 
ponents 

-V *kV = -V*Vk - k(V-V) . (2) 

The first term on the right represents horizontal advection of 
kinetic energy, while the second signifies flux due to horizontal 
divergence. Table 4 shows that both components produce an energy 
sink above 400 mb; however, outflow due to velocity divergence is 
the dominant process. 

Values of vertical flux divergence (Table 4) indicate an 
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export of energy in the lower and middle layers with significant 
import above 300 mb. This is consistent with the widespread 
ascending motion that characterizes much of the fifth SESAME 
period. The small total exists because vertical motions at 150 mb 
were adjusted to zero, whereas surface values were not. 

Current results now will be compared with those of Tsui and 
Kung (1977), hereafter denoted as T 8 K, Direct comparison with 
their study is possible since the resolution of their NSSL input 
data is very similar to that of the present investigation (meso 
8-scale). To the authors' knowledge, the energetics of no other 
periods have been investigated at such a fine resolution. (Mclnnis 
and Rung (1972) and Kung and Tsui ( 1975) considered the same days 
as T & K.) T & K computed energy budgets on nine days. These were 
categorized according to the particular type of storm conditions 
that occurred, i.e., convective, non-convective, or frontal. Four 
days were labeled as convective because active storms were identi- 
fied in the network along with strong, upward vertical motions. 
Three days were classified as non-convective since storms were 
either weak or failed to develop within the network area. 
Finally, two days were characterized by their association with 
cold fronts, either by the presence of a prefrontal squall line or 
by frontal passage through the analysis region. The development 
of convective activity in both cases, however, was limited. Table 
5 presents composite kinetic energy budgets for T 8 K’s categories 
along with current results. 
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Table 5. Mean kinetic energy budget in various types of storm situations (Tsui and Rung, 1977). 
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T & K’s convective category was the most energetically 
active. Vertical totals of generation on individual days ranged 
from 73.00 to 110.14 W m" 2 s The four day average was 89.45 IV m” 2 
(Table 5}. Dissipative losses tended to balance generation; verti- 

cal totals ranged from -14.08 to -200.44 IV ta t for a mean value 

-2 

of -95.99 W m , A significant amount of kinetic energy was tran- 
sported into the upper troposphere. Horizontal inflow existed in 
the middle and lower layers, whereas weaker values of flux diver- 
gence occurred in the 400-100 irb layer. 

Although a cold front was located in the mesonetwork during 
most of the AVE-SESAME V period, current energy results are very 

similar to those of T & K’s convective category (Table 5). The 

-2 

current vertical total for generation (83.57 W m ) agrees closely 

-2 

with their average value (89.45 W m ) . Ch the other hand, energy 
dissipation during AVE-SESAME V (-41.62 W m ) is approximately 
half that obtained by T & K (—95.99 W m” 2 ), but both represent 
significant losses either by traditional frictional processes or 
by transfers of energy to unresolvable scales of motion. Vertical 
profiles of these terms are similar for both studies, i.e,, large 
values near jet stream level, smaller values in the mid- 
troposphere, and secondary maxima near the surface. The current 
level of kinetic energy content (11.99 x 10 J m ) is approxi- 
mately half that of T & K (21.02 x 10^ J m” 2 ) , Although magni- 
tudes of local change are comparable, the AVESESAME V region 
experiences decreasing energy content, opposite that of T & K. 
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Finally, horizontal flux is a major sink for the current study, 
whereas it provides an important source for T & K. These con- 
trasts are probably attributable to differing locations of wind 
maxima between the two studies. 

Composite energetics of T & K’s non-convective category 
differs significantly from that of the convective periods and the 
current investigation (Table 5). The greatest contrasts are in 
the non-convective processes of cross-contour destruction and 
positive dissipation (energy supply from unresolvable scales of 
motion). Although energy processes for the frontal category agree 
in sign with those of the present study, frontal magnitudes are 
much weaker. This probably occurs because the convection during 
AVE-SESAME V was more intense than that during T & K’s frontal 
days. 


Since data were simultaneously collected at both the synoptic 
and meso 3 -scales during AVE-SESAME V, ve have a unique opportun- 
ity to compare the energetics of the same analysis region derived 
from different scales of input. Previous mesoscale energy inves- 
tigations (Mclnnis and Kung, 1972; Kung and Tsui, 1975; Tsui and 
Kung, 1977) did not have the corresponding 3 h synoptic- scale 
data; therefore, to the authors’ knowledge, this type of com- 
parison has never been previously reported in the literature. 
Using methods similar to those of Fuelberg and Jedlovec (1982), 
the synoptic- scale energy budget was evaluated at individual grid 
points for each observation time. Next, grid point values were 
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averaged within a 6.4 x 10 km area over Oklahoma which agrees 
closely in both size and location with the meso 3-scale region in 
Fig, 1. The resulting composite area-time averaged budget is 
given i n Table 6. 

There are many similarities between results from the two 
resolutions. Although energy processes dominating the meso 3 -scale 
storm environment usually characterize the larger scale flow, mag- 
nitudes at the finer resolution often are considerably greater. 
Closest agreements are found in energy content and its local 

derivative. The total for synoptic- scale kinetic energy content 
5 -2 

(13.20 x l(r J m ) compares favorably with the meso B-scale value 

c p 

(11.99 x ltr J m ) , and v ertical distributions show corresponding 
similarities. In addition, each area experiences a local decrease 
in energy during the 24 h period. 

Cross-contour generation (Table 6) is much larger at the meso 
3-scale (83.57 W m -2 versus 30.62 W m -2 ), probably due to the 
enhanced ability of the denser network to detect the strongly 
ageostrophic flow of the storm environment. Magnitudes of 

synoptic- scale dissipation are more comparable to those from the 
special sites. Horizontal transport at the synoptic scale is near 
zero. As a result, most of the energy generated within the region 
is transferred to unresolvable scales of motion. Ch the other 
hand, strong horizontal outflow is detected at the meso 3-scale, 
particularly in the 400-150 irb layer. Upward transport also is 
considerably stronger at the meso (3-scale. The greater magnitudes 
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of both flux terms at the smaller scale are due to at least three 


factors, First, since velocity divergence is highly scale depen- 
dent, horizontal transports at the finer resolution will be 
enhanced. This factor also produces greater vertical motions at 
the storm scale. Second, the special network better resolves 
small scale wind phenomena which are accompanied by strong hor- 
izontal gradients of kinetic energy, Finally, area-averaged 

energy content is greater at the meso g-scale. Thus, it appears 
that contrasts in results are mainly attributable to differences 
in resolution between the two data networks. 

Table 7 presents results from several previous studies of 
convective environments derived from synoptic- scale input data. 
Values from T & K's convective category and the current investiga- 
tion are included for easy reference. The energetics in the 

c o 

vicinities of two squall lines (areas approximately 11 x 10 km ) 
(Fuelberg and Scoggins, 1978) are most similar to current values, 
Generation by cross-contour flow is the greatest source in both 
studies while horizontal outflow and dissipation are major sinks. 
Investigations by Robertson and Smith (1980) and Fuelberg and 
Jedlovec (1982) considered much larger regions of the nation, and 
their results differ considerably from those of AVE-SESAME V, 
Judging from the wide variation in values from the larger areas, 
it appears that severe storm outbreaks can occur among many dif- 
ferent types of energy balances. This is attributable to the 
varying types of flow patterns associated with intense convection. 
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Table 7. Vertically integrated area-averaged kinetic energy budgets f 0 r various case stales 
Values of K are 1CT J m while others are W in . 



Synoptic 3ubrtu<;ioo 



Ch the other hand, results from AVE-SESAME V, T & K, and Fuelberg 
and Scoggins (1978) imply that the near storm environment may 
often be characterized by strong generation and negative dissipa- 
tion, regardless of whether synoptic or mesoscale data are 
employed . 

Previous studies have suggested that values of energy budget 
terms undergo major time variations that are related to the pres- 
ence and strength of nearby storms (Mclnnis and Kung, 1972; Kung 
and Tsui, 1975; Tsui and Kung, 1977; Fuelberg and Scoggins, 1978). 
Since intense convection occurred within the AVE-SESAME V network 
between 2000-0200 GMT, it is informative to examine current energy 
variability. Time series of surface- 150 mb totals and pressure- 
time cross sections are presented for this purpose. 

The time series (Fig. 10) shows that the meso $-scale region 

(Fig. 1) experiences a decrease of kinetic energy content through 

most of the 24 h period. The greatest drop occurs during the 

times of strongest convective activity (2000-0200 GMT) when values 

5 -2 

decrease 4.43 x 10 J m . The pressure-time cross section of the 
local change term (Fig. 11) indicates that most of this decrease 
occurs above 300 mb. Interestingly, however, a center of increases 
is located near 400 mb. This vertical variation will be discussed 
in greater detail in Section 4c. It agrees with findings of Vin- 
cent and Carney (1982) who observed decreases in jet level winds 
accompanied by the formation of a mid-level wind maximum during 
the Red River Valley tornado outbreak. T & K also observed a 
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Fig. 10. Time series of meso 8— scale energy budget terms 

integrated between the surface and 150 mb. Three 
hour time intervals are indicated on the horizontal 
axis. 
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Fig. 11. 


Pressure-t: 
locaj chanc 


W m /50 n± 



! cross section of area-averaged 
of kinetic energy. Units are 



buildup of mid-level energy prior to peak storm activity, but con- 
trary to current results, energy increases also were observed 
above 400 mb. The increases near 0800 GMT (Figs. 10-11) are due 
to jet intrusion into the network from the south. This feature 
will be described further in Section 4b. 

Vertical totals of cross-contour generation (Fig. 10) show 

major variations during the 24 h period. Values range from -11.6 

-2 -2 

W m at the initial observation time to a maximum of 208.98 W m 

at 2130 GMT. By 0500 CMT magnitudes drop to near pre- storm condi- 
tions. The pressure-time cross section (Fig. 12) shows that larg- 
est variability occurs above 400 mb and is centered about the time 
of greatest storm activity (2130 GMT). This extreme amount of 
upper level generation is attributable to strongly divergent flow 
being superimposed on a mesoridge (Fig. 6). Another interesting 
feature is the area of weak destruction that exists in the middle 
troposphere near 2130 GMT. T & K found significant energy genera- 
tion prior to and after thunderstorm passages, whereas cross- 
contour destruction existed within virtually the entire column at 
peak intensity. They suggested that generation surrounded the 
intense storm area while destruction occurred within the area. 
For the present SESAME case, however, destruction is limited to 
the middle layers only and is much weaker. Thus, the findings of 
T & K do not appear to represent all mesog-scale storm environ- 
ments. 

Figures 10 and 13 show that dissipation also undergoes major 
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fluctuations during the 24 h period. Greatest negative values 
occur in the upper troposphere centered about 2130 GMT. The 
secondary maximum at 0800 GMT is associated with the jet intru- 
sion. Patterns and magnitudes are quite similar to these of gen- 
eration (Fig. 12 ) r i.e, f transfers to unresolvable scales (nega- 
tive dissipation) correspond with cross-contour generation and 
vice versa. A similar agreement was observed by T & K Large 
negative dissipation suggests that much of the kinetic energy gen- 
erated at the meso g-soale by thunderstorm-induced environmental 
modification is transferred to scales of motion that are not ade- 
quately resolved by the special network. This transfer can either 
be to subgrid scales or to larger features which are poorly dep- 
icted by so small a data area. Finally, one should note the area 
of weak positive values in the middle troposphere near 2130 GMT. 
T & K observed a similar, although deeper and more intense, 
feature during their convective cases, and it has been described 
during synoptic- scale studies of convective environments as well 
(Vincent and Schlatter, 1979). 

Horizontal flux divergence provides a major energy sink (Fig. 
10) during the times of intense convection. The pressur — time 

cross section (Fig. 14) indicates that this process is greatest in 

—2 

the upper troposphere where values exceed -40 W m /50 mb. Con- 
sistent with the intense upward motion between 2000-2300 GMT, the 
cross section of vertical flux (Fig. 15) shows large upward tran- 
sport of energy. 
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Fig. 14. Pressure-time cross section of area-averaged 
horizontal flux divergence of kinetic energy 
Units are 10 * W m - 2/50 mb. 
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Fig. 15. Pressure-time cross section of area-averaged vertical 
flux divergence of kinetic energy. Units are 
W m~^/50 mb. 



Results of the area-averaged calculations demonstrate that 
the storm scale environment contains intense energy conversions 
and transports that are maximized at the time of strongest convec- 
tion. Table 8 quantifies conditions at 2130 GMT. As noted ear- 
lier, the environment is characterized by intense generation and 
negative dissipation in the upper troposphere. Horizontal outflow 
is a major upper level energy sink. In addition, present findings 
indicate that the energy balance is highly variable and fluctuates 
most rapidly during the times of storm activity. One should note, 
however, that several of the energy parameters began to change 
before the onset of the thunderstorms at 1700 GMT. Thus, some 
changes may be associated with atmospheric processes that helped 
create an environment suitable for the convection that followed. 
At additional factor may be diurnal variations in some energy 
terms (Kung, 1967). 

Finally, Fig. 16 shows time fluctuations of surface-150 rrb 
values over Oklahoma that were derived from the concurrent 
synoptic- scale data network. Except for horizontal flux, patterns 
of the terms are similar to those, derived from mesoscale data 
(Fig. 10); however, magnitudes from the NWS network are much 
smaller than those described earlier. The synoptic- scale network 
failed to capture the strong export of energy from the storm 
environment. In addition, energy variations associated with the 
jet intrusion at 0800 GMT are not indicated. 
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Table 8. Area-averaged energetics of the meso g^scale energy analysis region at 2130 GMT 
20 May 1979. Values of K are lCr J m while others are W m~ . 
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Fig. 16. Time series of synoptic-scale energy budget terms 
integrated between the surface and 150 nb for the 
energy analysis region (Fig. 1 ). Three hour time 
intervals are indicated on the horizontal axis. 
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b. Spatial distribution of energy parameters 


Horizontal fields of energy terms are presented to describe 
relationships between environmental wind patterns and various 
weather phenomena. Integrated kinetic energy content will be 
shown for all ten observation times, while the budget terms will 
be presented for two periods of interest. The first period 
involves energy variations associated with the convection (1700- 
0200 GMT), while the second concerns jet intrusion from the south 
near 0800 GMT. A 1 1 diagrams are for the 400-150 rrb layer since 
this is the region of most intense energetics. 

At 1100 GMT, a mesog -scale energy maximum is observed over 
southwestern Oklahoma (Pig, 17) . It moves northeastward and by 
1400 GMT is located just north of Hinton (HNT), At 300 mb (not 
shown), winds at Hinton are 48 m s compared to 36 m s for 
Clinton Sherman (CSM), situated only 80 km upwind. A significant 
new'feature at 1700 GMT is the energy decrease in the southeastern 
corner of the mesonetwork, It is also indicated by the local 
change term for 1700 GMT (Fig. 18) which, as one should recall, is 
obtained from centered differences. 

Energy content changes dramatically over northwestern 
Oklahoma between 1700-2000 GMT, coinciding with the period of 
rapid storm development in the southern portions of the network 
(Figs. 17-18). Values near Seiling (SSL) increase approximately 
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Fig. 18. Local change of kinetic energy for the 400—150 
mb layer. Units are W 
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Fig. 18. (continued) 
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40$, producing a localized energy maximum. Greatest values then 
decrease slightly by 2130 GMT, and by 2300 GMT, the maximum virtu- 
ally disappears. In contrast, energy content continues to 
decrease steadily in the southeast. Value, diminish nearly 50$ 
during the 6 h period ending at 2300 GMT. Energy content 
decreases significantly over most of the mesonetwork by 0200 CMT 
as the convection weakens. 

The timing and locations of the wind/energy variations sug- 
gest a relationship with the convective outbreak. Several previ- 
ous studies have noted the formation of upper level wind maxima on 
the poleward side of convective complexes (e.g., Fritsch and Mad- 
dox, 1981a; Fuelberg and Jedlovec, 1982; Fuelberg and Scoggins, 
1978). The current wind maximum in the northwestern corner of the 
data network also may have been convectively induced. Since 
mesoscale modelling studies by Fritsch and Maddox (1981b) revealed 
that a distinct speed maximum had formed after only one hour of 
simulated convection, the timing of the present case seems quite 
reasonable. A comparison with energy fields from the surrounding 
synoptic- scale network at 2000 GMT (Fig. 19) shows that the meso 
B-scale maximum over northwest Oklahoma may be the southern fringe 
of a larger jet feature located farther north. 

The wind/energy fluctuations in the southeastern part of the 
network also may be convectively related. However, a somewhat 
puzzling aspect of AVESESAME V is that these upper level winds 
began to decrease between 1400-1700 GMT, prior to the onset of 
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convection. This implies that storm- environment interactions were 
not the sole factors producing the wind variability. In fact, the 
energy diagram from the synoptic- scale network at 0200 GMT (Fig. 
19) reveals that the minimum „ T rsr southeastern Oklahoma is part of 
a much larger feature centered over Arkansas. Nevertheless, those 
changes occurring after 1700 GMT may have been convectively influ- 
enced to some extent. In addition, Section 4c will describe a 
vertical wind couplet in this region that appears to be storm 
induced . 

Patterns of the remaining terms now will be described just 
prior to and during the convective outbreak. Pronounced changes 
in cross-contour generation are observed between 1700-0200 GMT 
(Fig. 20 ). At 1700 GMT, maximum generation occurs just east of 
the Texas panhandle. Values increase significantly by 2130 GMT 
when greatest generation is found northeast of the convection. 
Strongly divergent and ageostrophic flow with large crossing 
angles (see Fig. 6) produce generation in excess of 300 W m . 
Values decrease dramatically by 0200 GMT as the diffluence and 
mesoridging diminish in response to the weakening convective 
activity. 

Velocity divergence at 1700 GMT (not shown) is relatively 
weak; therefore, patterns of horizontal flux divergence at this 
time (Fig. 21) mainly reflect advection of kinetic energy. A com- 
parison with the corresponding diagram of energy content (Fig. 17) 
shows that inflow is found downwind (northeast) of energy maxima 
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and vice versa. Flux values strengthen to as much as -300 W a 
by 2130 GMT when strong upper level velocity divergence and hor- 
izontal advection produce energy export over most of the region. 
The only exception is the northcentral area of energy import which 
is attributable to velocity Convergence. Values diminish by 0200 
GMT. At 2130 GMT, upward vertical motion reaches -60 yb s at 
500 rrb in the southeastern edge of the region (not shown). This 
produces large import of energy from the lower levels (Fig. 22). 
The much weaker vertical export in the northcentral region may be 
due to compensating downward motion outside the near storm 
environment ■ 

Finally, fields of dissipation (Fig. 23) show that much of 
the energy generated within the storm environment is transferred 
to non-resolvable scales of motion, thereby providing a prominent 
energy sink to the meso B-scale. At 2130 GMT, however, the south- 
ern portion of the network containing the storms has the smallest 
negative values and, in fact, exhibits a small area of positive 
dissipation. The primary area of negative values is located 
farther north, near the region of greatest cross-contour genera- 
tion (Fig. 20) . 

A second interesting feature during AVESESAME V is the 
intrusion of a wind maximum into the special network from 'the 
south during the final third of the experiment period (Figs. 17- 
18). The perturbation is first evident at 0500 GMT, and maximum 
energy content is reached at 0800 GMT when levels had increased 
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Fig. 22 . (continued) 
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Fig. 23. Dissipation of kinetic energy for the 400-150 mb 
layer. Units are 10^ W 
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Fig. 23 . (continued) 
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approximately 50$ from those 6 h earlier. By 1100 GMT 21 May, it 
had weakened and/or moved out of the special network. Mechanisms 
producing the wind maximum are unknown; however, it may be related 
to thunderstorm activity in northcentral Texas. The feature w'-us 
not detected by the NWS stations. Spatial maps for 0800 GMT show 
the energy balance for the feature, Horizontal flux convergence 
exceeding 400 W m is its greatest energy source (Fig. 21). Cal- 
culation of the two flux components (see (2)) indicates that 
energy import primarily results from advection, and not velocity 
convergence, Another important energy source to the southcentral 
region is cross-contour generation (Fig. 20) which reaches 300 W 
m . A center of destruction is located downwind of the energy 
maximum, apparently due to supergradient flow in advance of the 
feature . 

Dissipation due to unresolvable scales of motion that reaches 

—2 

-700 W m is the greatest energy sink to the area of the wind 
maximum (Fig. 23). Farther north of the feature, however, a 
region of weaker positive dissipation is located near the area of 
supergradient flow (negative generation). Previous synoptic- scale 
studies also have noted destruction of kinetic energy and positive 
dissipation in advance of propagating jet streaks (e,g, f Robertson 
and Smith, 1980). Finally, downward motion provides a small 
energy sink to the jet area (Fig. 22). 


j 


78 



c. Vertical variations in wind 


Horizontal maps of kinetic energy content revealed decreases 
in the southeastern portion of the network between 1400-2300 GMT 
(Fig. 17). Mc*e detailed diagnoses in this section indicate a 
complex vertical pattern that is superimposed. Cross sections of 
energy content along an axis from near Seiling (SEL) southeastward 
to near Elmore City (EMC) are presented in Fig. 24. The axis 
passes through both the energy maximum and minimum seen earlier 
(Fig. 17). Greatest energy content at 1700 GMT is located in the 
northwest, near 250 mb. Between 1700-2000 GMT, upper level values 
increase near Seiling and Hinton. Ch the other hand, significant 
energy decreases occur in the southeast, most notably in the BOO- 
200 mb layer. Energy content then decreases between Seiling and 
Hinton by 2300 GMT, but a vertical couplet develops in the 
southeast, where a minimum was observed in the 400-150 mb layer 
(Fig. 17). Notable energy increases occur near 400 mb while a 
significant energy decrease is seen near 200 mb. Although most 
pronounced over the southeastern part of the network, the feature 
was detected in the area-averaged local change term (Fig. 11), 
thereby indicating that it is characteristic of a rather large 
part of the AVE-SESAME V region. 

The time series of horizontal wind speed for Healdton (HEA) 
clearly shows the fluctuations (Fig. 25). Above 300 mb, speeds 
decrease slowly between 1100-2000 GMT, but then decrease nearly 
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50$ during the 3 h period between 2000-2300 GMT. Ch the other 
hand, values nearly double near 400 irb in the same 3 h period. As 
a result of the rapid fluctuations, the altitude of maximum wind 
speed drops approximately 200 mk. in only 3 h. By 0500 GMT, the 
wind patterns return to near original configurations. Similar 
variations near 2300 GMT also appear in the time series for Chic- 
kasha (CHK), Mountain View (MTV), KIW (TVY), Hinton (HNT), 
Wichita Falls (SPS) , Altus (LTS), and Fort Sill (FSI), although 
they are not as pronounced. These agreements suggest that the 
changes are valid and not merely attributable to data uncertain^ 
ties. 


The upper level wind variations appear to be related to the 
nearby convective activity. Similar fluctuations have been docu- 
mented during AVE-SESAME I which coincided with the Red River Val- 
ley tornado outbreak (Vincent and Carney, 1982; Vincent & 
1981). In that case, a mid-level jet formed near 500 nh during the 
thunderstorm activity. Upper level winds showed corresponding 
decreases. When the convection diminished, the mid-level jet also 
dissipated . 

To isolate energetics of the wind couplet from that of the 
remainder of the analysis region, grid point values of the energy 
budget were averaged over a fixed subvolume (1.5 x 10 4 km 2 area) 
in which the feature was most pronounced (Fig. 1). The pressure- 
time cross section of the local change term (Fig. 26) shows that 
the couplet occurs between 2000-2300 GMT. Energy decreases are 
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found in the 350-150 irb layer, whereas greatest increases occur 
between 550-350 mb. Table 9 presents integrated totals of the 
energy budget for both of these layers. 

Energy content within the 350-150 rrb layer decreases signifi- 
cantly during the 3 h period of intense storm activity. The ori- 
ginal value of 7.63 x lCT* J m" 2 at 2000 GMT drops nearly 50% to 
3.95 x 10 5 J m ~ 2 at 2300 GMT. Table 9 and Fig. 26 indicate the 
energy mechanisms producing the vertical couplet. With the excep- 
tion of large values near 0800 GMT which are due to the jet intru- 
sion described earlier, greatest magnitudes occur near 2130 GMT. 
This finding is similar to that observed for the entire analysis 
area (Fig. 10). At 2130 GMT, generation and vertical flux conver- 
gence provide sources to the 350-150 irb layer that are 136.50 W 
-2 —2 

m and 95.90 W m , respectively. However, energy sinks by hor- 

m 

—2 

izontal outflow (-130.45 W m ) and negative dissipation (-136.00 

—2 

W m ) override to produce an overall energy decrease. The rela- 
tive contributions of the four budget terms are similar at all 
three times of intense convection (Table 9). 

Contrary to findings of the topmost layer, energy content 
between 550-350 irb increases during the same 3 h period (Table 9 f 
Fig. 26). At 2000 and 2130 GMT, significant sinks by vertical 
flux divergence and weaker sinks by horizontal export are coun- 
teracted by sources due to generation and positive dissipation. 

It is noteworthy that at 2000 GMT, energy transfer from non- 

—2 

resolvable scales (33.15 W m ) is an even greater source than 
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that due to cross-contour flow (5.69 W m ). At 2130 GMT, posi- 
tive dissipation is only slightly smaller than generation. The 
energy balance changes somewhat by 2300 QvT when generation and 
horizontal import are in near balance with negative dissipation 
and vertical flux divergence. 

Results indicate that the complex vertical fluctuations of 
wind (energy) are due to differing balances between the various 
budget terms. The importance of positive dissipation is espe- 
cially significant. Since dissipation is computed as a residual, 
its uncertainty is somewhat greater than that of the other terms. 
Nonetheless, sensitivity analyses demonstrate that errors are not 
its dominant component (e.g., Appendix, Ward and Smith, 1976; 
Fuelberg and Jedlovec, 1982). Large magnitudes of dissipation at 
times of greatest convective activity ( 2000-2300 GMT) suggest that 
the near storm environment is highly influenced by motions which 
are not adequately resolved by the meso 0-scale data network. One 
possible explanation is that meso y- scale vertical motions within 
the deep convection are "mixing” the air, thereby producing a sink 
at the original jet level and a corresponding source at levels 
below. This effect would not be measured by the other budget 
terms which describe meso 8 - scale mechanisms. In a somewhat simi- 
lar manner, Vincent and Schlatter (1979) noted that subgrid scale 
processes were an energy supply to the 500-300 mb synoptic- scale 
flow during times of intense convection. They speculated that 
potential energy on the cumulus scale was transferred to kinetic 
energy on the synoptic scale, 
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5. SUMMARY AND CONCLUSIONS 


A meso 8- scale kinetic energy analysis was conducted for a 
period of intense thunderstorm activity that occurred over central 
Oklahoma during AVE-SESAME V (20-21 1V% 1979). Rawinsonde data at 
75 km spacings at either 3 or 1.5 h intervals were used to examine 
energy processes of the mesoscale storm environment prior to, dur- 
ing, and after peak convection. The goal was to better understand 
how areas of storms modify their surrounding larger scale environ- 
ment. Unlike previous meso 8- scale energy studies which 

emphasized a collection of days, the current investigation was a 
detailed analysis of a single 24 h period. 

Area-averaged energetics indicated that the meso 3 -scale 
storm environment was characterized by generation of kinetic 
energy by cross-contour flow, horizontal flux divergence of 
energy, upward transport by ascending motion, and transfer of 
energy from resolvable to non-resolvable scales of motion (nega- 
tive dissipation). Vertical profiles showed that processes were 
greatest near jet stream level. 

Current results were compared with previous studies by Tsui 
and Kung (1977) (T & K) who computed energy budgets for nine days 
using meso 8- scale NSSL data. Findings of AVESESAME V were very 
similar to those from T & K's composite convective group. Both 
were characterized by strong generation and significant energy 
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losses by negative dissipation. Although horizontal flux diver- 
gence was a major energy sink during AVE-SESAME V, this was not 
the case in T & K’s investigation. Differing locations of wind 
maxima between the two studies may have contributed to this con- 
trast. 

A unique feature of AVE-S1SAME V is that rawinsonde data were 
simultaneously collected at both the synoptic and meso @ -scales. 
Therefore, a comparison of the kinetic energy balance of the same 
analysis region based upon different scales of input data received 
special attention. To the authors' knowledge, such a comparison 
has never previously been reported in the literature. The ener- 
getics of the two resolutions was similar in that generation and 
dissipation dominated both the mesoscale storm environment and the 
synoptic- scale flow. However, magnitudes were larger at the finer 
resolution. The coarser network failed to detect the enhanced 
horizontal flux divergence within the storm environment. 

The meso 3 -scale energy budget underwent major time fluctua- 
tions. Pressure-time cross sections showed that maximum variabil- 
ity occurred above 400 mb, centered about the time of greatest 
storm activity. Strong upper level generation during the convec- 
tive outbreak was attributable to storm-induced divergent flow 
being superimposed on a mesoridge. Magnitudes and variability of 
horizontal flux divergence and dissipation were very similar to 
those of generation, thereby indicating that much of the kinetic 
energy generated in the mesoscale storm environment was exported 
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outside the region and transferred to non-resolvable scales of 
motion. Except for horizontal flux, time variability of energy 
parameters derived from synoptic- scale data was similar to that 
from the special network. 

Horizontal fields of energy content for the 400-150 rrb layer 
were presented to describe flow patterns of the meso 6-scale storm 
environment. Energy content was seen to change dramatically dur- 
ing the period of intense convection. Winds strengthed abruptly 
north of the convective area as the severe thunderstorms 
developed. Corresponding values of energy content increased 
approximately 40$ in only 3 h. These findings agree with those of 
jet development poleward of Mesoscale Convective Complexes. In 
contrast, energy content southeast of the storms decreased slowly 
prior to the convective outbreak, but then diminished as much as 
50% once the convection developed. In addition, the horizontal 
maps showed the intrusion of a wind maximum into the region from 
the south during the latter part of the experiment period. Energy 
levels in that area increased approximately 50$ during a 6 h 
period. 

Horizontal maps of the energy budget parameters showed good 
agreement with observed weather features. During the times of 
intense convection, strong cross-contour flow produced large gen- 
eration north of the storm area. Fields of dissipation indicated 
pronounced transfers of energy to non-resolvable scales over most 
of the region. Negative centers were co-located with areas of 
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maximum generation. Interestingly, however, at the time of peak 
storm intensity, a small area of positive dissipation existed near 
the convection. Strong horizontal flux divergence was an energy 
sink over most of the region, whereas vertical flux convergence 
was a major source. The area of jet intrusion toward the end of 
the period was characterized by horizontal flux convergence, 
cross-contour generation, and negative dissipation, 

Finally, detailed diagnoses at the time of peak convection 
revealed a complex vertical pattern of wind variations above 600 
rrb in the southeastern part of the network. The feature was super- 
imposed on generally decreasing winds that began prior to the 
storm outbreak. Speeds near 200 nb decreased approximately 50$ 
during the 3 h period coinciding with the most active storms. Ch 
the other hand, winds near 400 rrb nearly doubled during the same 3 
h period. As a result, the level of maximum speeds dropped 
approximately 200 mb. To explain this feature, a kinetic energy 
analysis was conducted for a subvolume of the total analysis 
region to isolate the speed variations. Results indicated that 
the vertical fluctuations of wind (energy) were due to differing 
balances among the various energy budget terms. Large magnitudes 
of dissipation suggested that the mesoscale storm environment was 
highly influenced by motions that were not adequately resolved by 
the meso £ -scale network. Vertical variations in dissipation con- 
tributed greatly to the vertical variation of wind. 


Energy/wind relations during AVE-SESAME V appear consistent 



with those of previous storm- environment studies. However, AVE- 
SESAME V is a singularly well documented case. W hypothesize 
that the thunderstorms produced many of the observed energy varia- 
tions. Scale interactions appear extremely complex in the near 
storm environment. Unfortunately, our understandihg of sub- 
synoptic phenomena, as well as scale interactions, is very lim- 
ited. Thus, many more theoretical, diagnostic, and numerical 
investigations are needed to explain the observed phenomena. As 
our overall knowledge of the smaller scales improves, greater 
forecasting skill will result. 


92 


a 



APPENDIX 


Sensitivity Analysis 

The purpose of the sensitivity analysis is to quantify the 
effects of random errors in the rawinsonde data on the kinetic 
energy budget parameters. ,0ur procedure was to add computer gen- 
erated random perturbations that simulated actual rawinsonde 
errors to the original time adjusted 25 mb data at individual sta- 
tions. Then, energy budgets were calculated from the perturbed 
data and compared with those derived from the original data. The 
sensitivity is represented by mean absolute differences and mean 
correlation coefficients between the original and perturbed budg- 
ets. This technique is similar to that used previously by Ward 
and Smith ( 1976), Robertson and Smith (1980), and Fuelberg and 
Jedlovec (1982). 

Perturbations introduced into the original data were normally 
distributed about zero with standard deviations as a function of 
pressure (Table 10). Deviations for levels not shown can be 
obtained by linear interpolation. The values are similar to those 
suggested by Fuelberg and Jedlovec (1982). Because of the care 
taken in processing and checking the data, perturbations were res- 
tricted to two standard deviations from zero which included about 
95 % of the possible values. 

Since maximum convective activity and strongest energetics 
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Table 10. Standard deviations of perturbations normally 
distributed about zero. 


Pressure 

level 

(mb) 

Wind 

direction 

(deg) 

Wind 
speed 
(m s -1 ) 

Height 

(m) 

150 

10.0 

4.3 

32.0 

200 

9.1 

4.0 

28.0 

300 

7.0 

3.7 

21.0 

500 

4.6 

3.0 

11.0 

700 

3.5 

2.0 

8.0 

900 

2.0 

1.0 

4.1 


i 
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occurred near 2130 GMT 20 May, this observation time was chosen 
for the analysis. Ten versions of the data, each with a different 
set of perturbations, were made for 2130 GMT. Analytical pro- 
cedures (Barnes analysis, vertical motions, etc.) used with the 
perturbed data were identical to those used with the original data 
(see Sections 2c and 2d). By comparing the original and perturbed 
data grids, it was found that approximately 35 % of the magnitudes 
of height perturbations introduced at the station locations (Table 
10) were removed by the Barnes analysis scheme. Since winds were 
vertically filtered before the objective analysis, approximately 
55 % of their introduced error was removed. Thus, our computa- 
tional procedures have the desired ability to reduce the inherent 
error within the data set. 

Original budget values along with mean absolute differences 
and correlation coefficients between original and perturbed ver- 
sions are given in Table 11 for the entire energy averaging region 
and in Table 12 for the limited area in the southeast (seeFig.1). 
Results show that the dissipation and generation terms are most 
sensitive to random error. The greater uncertainty in dissipation 
is expected since it is computed as a residual from the other 
budget terms (see (1)). Mclnnis and Kung ( 1972) have commented on 
the difficulties in accurately calculating the generation term. 
Ch the other hand, terms representing energy content and vertical 
and horizontal flux divergence appear to be the most reliable. 
With the exceptions of dissipation for the 700-400 rrb layer over 
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Table 11 Area-averaged kinetic energy budget for 2130 GMT 20 May 1979. Values in parentheses a rt» 
mean absolute differences between the original and ten perturbed budgets. Values in 
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Vertical [0.99] Q.97] [0.82] [0.95] [0.99] [0.88] 

Total (0.27) q.82) (25.46) (7.92) (0.00] (26.83) 



S3 cm 

° ‘(3 


T 3 

5 

!5 


cfl i-j 


05 
(U 
O 

53 ■ ■ 
a) cum 
ft o 

ft 
o 

o co 
•H 
S 3 

cfl 43 
" cJ 


ft 
<U 

rH 44 
Cfl 44 
O *H 

•rl Tfl 
44 -1 

ft 0 ) 0 ) 

cu c-» e 
0 ) o 


A 
CU & 
ft 


43 co cfl 64 
4 J 4 = 


Cfl CO ft 
4-1 O 
0 ) 44 


a (3 - 
fi Cfl 
co 41 o 4J 

o B cfl ft 

H ft CU 
o a) 43 

C ft 

8 cfl q 
°h 

<U Cfl *> 

■ r4 C0T-4 
CO <U l _ 

, . <D 3> S 

O 43 H 

a) <u 

ft ft 
cfl * cfl 
a) &. co 

r* 44 CO 

■ ■COM 
•H 00 <U 

T3 


I 


■§ 

co 


o 

0) 


4-4 

o 


to fl fl 


a) 43 

CO pi 
Hrffl 


co cfl a) cfl 
>>> rO & 


cfl ft 
r -4 3 

• 44 
O /-s >4 

g i-t <u a) 
4J & q 


00 

ft 


M M fi 
O 'H 4) 
44 pH 4-> 


0) 



>s< 

60 

>4 O 

ai cn 

rt H 

a) cm 

(0 

T3 44 43 
ai cci 4-i 
44 

Cfl 4-1 0 

m ai a) 

60 t— 4 a) 44 
<0 ft S 44 
y 3 4j ai 
pi o cu o 

H 04 CJ 


04 


co 

4-1 

(3 

CU 

•H 

CJ 

•H 


■s 

H 


O- 

<ro 

3 

co 


% 

fl> 


O 

o r-- m 
• -fl- 

\Q • • 

CO o co 

r-H 1 » *— < 

I 


O f *"“ 1 

a\ cr* «n 

ON 1 — 1 

m • * 

ON o oo 


-e- 

■t> 


lO « — * /*s 
<r r- on 
• on 

o • • 

CO O N 

I 4—4 |— 1 

I ^ 


O ' — 1 

m <r 
• co f* 
• • 
co o r^. 


■p 

co 

>? 

CO 




m i— • 

o m c 

• <T» O 
<t • 

CO O C 

| i— 


vO i— > / 
<TNC 
• ON r 
in * 
O c 


CU <-> 
ft A 
PJ 6 
co ^ 
co 


CU 

ft 

P4 


ft 

CU 

> 

cfl 

A 


O 

m 

T— 4 

I 

o 

m 

cn 


IS- l—l p— v 

<y> in o 

oo • 
mom 


rs p — 1 - 

m oo t-t 

. 05 m 
\Q • • 

vD O 

| l i v--' 


CM r-n/-s 
VO O lO 
• ON CO 

i— i • • 

I o CO 


00 r— i /-> 

vD CO CO 
. 00 CM 
ON * * 

CO O CM 


VO r — 1 

<r on m 
• on <r 

o • • 
o o 


on « — i 
in on co 

• O H 
CM • • 

O O 


o 

m 

co 

o 

m 

m 


97 


J 


the analysis area, and horizontal flux for the 550-350 irb layer 
over the limited area, percentages of mean absolute differences to 
original budget values are less than 35 % (not shown). At any 
particular level, one can have greater confidence in terms having 
large magnitudes than those near zero. 

To quantify the sensitivity of the horizontal energy pat- 
terns, linear correlation coefficients were calculated between 
grid points of the original field and each of the ten perturbed 
runs. Results of these calculations also are given in Tables 11 
and 12 . With the exceptions of the generation and dissipation 
terms, mean coefficients for the ten runs are never less than 0.90 
for either area. For the limited region, the dissipation term of 
the 350-150 irb layer has the lowest mean correlation (0.77); 
whereas for the entire region, the generation term has the lowest 
correlation (0.80) in the 700-400 nb layer. One should note that 
the dissipation term has a high correlation ( 0 . 95 ) in the 550-350 
irb layer of the limited area. This is the layer where subgrid- 
scale energy transfers were a significant energy source. 

Horizontal fields having the lowest correlations of the ten 
perturbed runs are given in Fig. 27 for generation, horizontal 
flux, and dissipation of energy, respectively. Original fields of 
these terms are given in Figs. 20 , 21 , and 23 . Generally, pat- 
terns of the worst perturbed fields show good agreements with 
their corresponding originals. Patterns are nearly identical for 
horizontal flux divergence, but somewhat greater discrepancies can 
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be seen for the other terms. Although data errors sometimes pro- 
duces considerable variations at individual grid points, horizon- 
tal patterns usually contain the same major features. 

In summary, results of the sensitivity analysis show that 
random errors in rawinsonde data generally should not affect the 
interpretations of the current study. In addition, good space and 
time continuities of the original patterns give additional confi- 
dence in final results. 
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